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Using numerical simulations we examine the velocity fluctuations of a probe particle driven with 
constant force through a two-dimensional disordered assembly of disks which has a well defined 
jamming point J at a density of <f)j — 0.843. As cj> increases toward (j>j, the average velocity 
of the probe particle decreases and the velocity fluctuations show an increasingly intermittent or 
avalanchelike behavior. When the system is within a few percent of the jamming density, the velocity 
distributions are exponential, while when the system is less than a percent away from jamming, the 
velocity distributions have a non-exponential or power law character. The velocity power spectra 
exhibit a crossover from a Lorentzian form to a 1// shape near jamming. We extract a correlation 
length exponent v which is in good agreement with recent shear simulations. For <f> > <j>j, there 
is a critical threshold force F c that must be applied for the probe particle to move through the 
sample which increases with increasing (f>. The onset of the probe motion above <f)J occurs via a 
local yielding of the particles around the probe particle which we term a local shear banding effect. 

PACS numbers: 45.70.-n,45.70.Ht,05.40.Ca 
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I. INTRODUCTION 

When a disordered assembly of particles with short- 
range steric interactions becomes dense enough that the 
particles begin to overlap, a transition occurs from an 
unjammed state where the system easily flows or acts 
like a liquid to a jammed state where the system resists 
shear and shows a solid-like response. Liu and Nagel 
proposed a general phase diagram for disordered particle 
assemblies where a transition into a jammed state occurs 
as a function of several different variables including in- 
creasing density, decreasing external forcing, or decreas- 
ing temperature At 'Point J' on this phase diagram, 
jamming occurs for increasing density at zero stress and 
zero temperature |2|]. There is growing evidence that 
point J has the properties of a critical point with diverg- 
ing length scales on both sides of the jamming transition 
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has been well studied in the particular system of two- 
dimensional, finite range repulsive frictionless bidisperse 
disks with a radii ratio of 1 : 1.4 @, H [13, E d. The 
bidispersity is necessary to prevent the system from crys- 
tallizing. Simulations show that the jamming density 
« 0.84, and experimental realizations of this system 
also give a similar jamming density flol Il6j. 

One aspect of jamming that has been less well studied 
is the dynamical response as point J is approached. If 
point J is a critical point, it could have strong effects on 
the dynamical response of the system to global or local 
forcing and may give rise to intermittent or avalanche 
behavior. In the numerical work of Drocco et al. Q, the 
velocity fluctuations of a single probe particle driven with 
constant force showed signs of increasing intermittency 
and avalanche characteristics as (j) approached </>j. More 
recent experimental work in which a single probe particle 
was driven at constant force through a granular packing 
also indicates that the velocity fluctuations become in- 



creasingly intermittent as (f)j is approached [16j. These 
results suggest that proximity to point J will strongly af- 
fect the dynamics and may result in the appearance of 
critical properties. 

The velocity and force characteristics of a locally 
driven particle have already been used as a probe of the 
yielding and dynamical properties of a two-dimensional 
disordered colloid system, where a finite threshold force 
must be applied before the particle can be set in mo- 
tion [ij], [3. Experiments and simulations of three- 
dimensional disordered colloidal systems have shown that 
a threshold force for motion appears when the de nsity 
is high enough for a glass transition to occur plj l'2[| . 
while for jamming in two dimensions a finite threshold 
force is also expected to appear above the jamming tran- 
sition. Drag effects in granular systems have been studied 
by moving a large obstacle through the granular media 
pH |22| . Geng and Behringer analyzed the fluctuations 
of a probe particle held fixed in a two-dimensional sys- 
tem while the remaining grains were driven past [23| , and 
found that the force fluctuations have intermittent char- 
acteristics and are exponentially distributed. 

In this work we examine the velocity fluctuations of 
a probe disk driven with constant force through a two- 
dimensional bidisperse system of disks as the density of 
the system passes through point J. As <j)J is approached 
from below, the velocity fluctuations show increasingly 
intermittent properties characterized by large bursts of 
motion. The average probe particle velocity initially de- 
creases linearly with increasing (f>; however, within a few 
percent of the jamming density, the average velocity falls 
to zero in a nonlinear manner. The velocity fluctuations 
are exponentially distributed when <f> is within several 
percent of the jamming density, while within one percent 
of 4>j the velocity fluctuations have a non-exponential 
distribution which can be fit with a power law. The 
power spectra shows a Lorentzian behavior with al/w 2 
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behavior at high frequencies that is similar to the results 
found in the slow drag experiments of Ref. (23| ; however, 
very close to jamming the low frequency response ex- 
hibits al/w behavior. The local rearrangements of disks 
near the probe particle increase in spatial extent as the 
jamming transition is approached from below. 

For <j> < (f)j, the probe particle is always able to move 
through the sample for arbitrarily small driving forces. In 
contrast, for cf> > <f>j, there is a finite threshold force F c 
required to move the probe particle. Above the thresh- 
old force, progress of the probe particle is made possible 
by yielding of the surrounding particles that occurs via 
local plastic rearrangements of the disks near the moving 
probe particle. The appearance of a finite driving thresh- 
old in the jammed phase is similar to recent numerical 
studies of dense colloidal systems where a threshold force 
appeared in sufficiently dense samples [2(| ■ We find that 
the threshold force increases with increasing <j). The ve- 
locity fluctuations at drives just above the yielding tran- 
sition F c have a nonexponential distribution, while at 
higher drives the velocity fluctuations become exponen- 
tially distributed. When the probe particle moves, it 
creates local rearrangements or a local shear banding. 



II. SIMULATION 

We simulate a two-dimensional bidisperse system of N 
disks of size L x L with periodic boundary conditions 
in the x and y directions. The particles arc modeled 
as repulsively harmonically interacting disks of radii Ra 
or Rb with a size ratio Ra ■ Rb of 1.4 : 1. The ratio 
of A and B particles is 50 : 50. The initial disk pack- 
ing is obtained by placing the disks in non-overlapping 
locations and then shrinking all disks, adding a few ad- 
ditional disks, and reexpanding all disks while thermally 
agitating the disks. This process is repeated until the de- 
sired density is obtained. After the disk configuration is 
generated, we select a single disk and drive it with a con- 
stant force Fjj = FjjU. The motion of an individual disk 
i at position is governed by the following overdamped 
equation: 

^ = E k ( R °ff " M^M^ff - |r«|) + F l D 

Here 8 is the Heaviside step function, r,j = rj — r,-, 
and R e ff = Ri + Rj, where Ri and Rj are either Ra 
or Rb depending on the type of each disk. We take the 
spring constant k — 200. The driving term = Fdx is 
applied to only one disk and we can extract the instanta- 
neous velocities of all disks. In this work, L — 60 and the 
jamming density is <fij w 0.843, corresponding to approx- 
imately N ~ 2600 disks. The ability to experimentally 
realize a system where the probe particle velocity can 
be examined under a constant driving force has recently 
been demonstrated [16j |. 
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FIG. 1: Images of the system at different densities where 
the probe particle (red) is driven 2/5 of the way through a 
disordered sample of bidisperse disks. The smaller disks are 
light colored (green) and the larger disks are dark colored 
(blue). Black lines indicate the trajectories of the disks that 
have moved, (a) <j> — 0.13, (b) (f> = 0.323, (c) <f> = 0.55, and 
(d) = 0.81. 



III. VELOCITY FLUCTUATIONS BELOW 
JAMMING 

In Fig. 1 we show snapshots of the system at differ- 
ent densities of </> = 0.13, 0.323, 0.55, and 0.81. The 
probe particle is highlighted and the trajectories of all 
grains are indicated as the probe moves 2/5 of the way 
through the sample. The larger disks are dark colored, 
the smaller disks are light colored, and the overall disk 
arrangement is disordered. At each of the densities pre- 
sented, the probe particle leaves an empty space behind 
it after pushing the other disks out of the way, as is most 
evident in Fig. 1(c). The appearance of an empty trail 
behind the probe particle generally occurs for any density 
below (f>j since there is no pressure in the system to force 
the surrounding particles back into the empty space. 

To construct velocity histograms for different densities 
at fixed driving force, we perform a series of simulations 
at increasing values of <f> and measure the instantaneous 
time trace of the probe particle velocity in the direction 
of the drive, V(t). We then average the velocities over 
the entire trace to obtain (V). In Fig. 2(a) we plot (V) 
vs cj) for a system with Fo = 0.025. Here (V) monotoni- 
cally decreases with </> and goes to zero at <pj — 0.843, a 
jamming density which is in agreement with other stud- 
ies [H, Fig. 2(a) also shows that slope of the de- 
crease in (V) with increasing <fi becomes shallower above 
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FIG. 2: (a) The average velocity (V) of the probe particle vs 
4> for Fd — 0.025. At low densities, (V) varies linearly with 
4>, and a small change in the slope appears near <f> = 0.5. (V) 
goes to zero at the jamming density <pj — 0.843, indicated 
by the letter J. (b) A detail of (V) vs cf> from (a) showing 
that {V) does not start to fall rapidly until 4> is within three 
percent of the jamming density. 

4> = 0.5, and that (V) changes from a roughly linear de- 
crease with increasing to a much more rapid nonlinear 
decrease for <j> > 0.82. When <p < 0.5, the probe par- 
ticle motion alternates between free motion unimpeded 
by any disks and occasional collisions during which the 
probe particle pushes one additional disk along with it. 
We show later that this appears as a clear feature in 
the velocity distributions. For (f> > 0.5, the probe par- 
ticle spends most of its time pushing one or more parti- 
cles but the number of dragged particles remains small. 
The distance over which the probe particle drags indi- 
vidual background disks varies widely due to interactions 
among the dragged disks which create random fluctua- 
tions for the dragged disk closest to the probe particle. 
For <fi > 0.82, the induced motions in the background 
packing become increasingly collective, and even disks 
which are far away from the probe particle can partici- 
pate in the motion. In Fig. 2(b) we plot (V) versus <j> 
close to jamming. The velocity decreases in a nonlin- 
ear fashion within three percent of the jamming density. 
Within less than one percent of (f>j, the velocity curve 
is still nonlinear; however, the concavity of the slope of 
(V) versus (f> changes sign. In a later section we discuss 
how the velocity scales near jamming and how it can be 
related to a diverging correlation length. 

We next examine in detail the velocity distributions 
obtained from time series V(t) taken at fixed values of <f). 
In Fig. 3(a) we plot V(t) for a sample with <j> = 0.323 and 
Fd = 0.025, the same system shown in Fig. 1(b). This 
single time trace is taken during the period required for 
the particle to move slightly more than halfway through 
the sample. Figure 3(a) shows that the probe particle of- 
ten moves at the maximum possible velocity V = 0.025. 
The plateaus at V = 0.025 occur when the probe par- 
ticle is not in contact with any other disks. Collisions 
with other disks reduce the velocity of the probe particle 
and in some cases the probe pushes one or more back- 
ground particles for a period of time, creating character- 
istic drops in V(t) of 50% or more as seen in Fig. 3(a). In 
Fig. 3(b) we plot the distribution of the velocities P(V) 
for the system in Fig. 3(a). The histogram is obtained 




FIG. 3: (a) V(t), the time trace of the probe particle ve- 
locity, at (f> — 0.323 for Fd = 0.025. Between collisions the 
probe particle moves at a velocity of V — 0.025 due to the 
overdamped dynamics, (b) The histogram P(V) of the ve- 
locities for (f> — 0.323 for ten different initial conditions at 
F D = 0.025. The peak at V = 0.025 corresponds to the 
collision-free portion of the particle motion, while the sec- 
ond peak near V = 0.0125 corresponds to the times that the 
probe particle is pushing one extra particle which doubles 
the drag. All of the P(V) in this work are normalized such 
that the probabilities sum to one. (c) V(t) at (j> — 0.55 and 
Fd — 0.025. (d) P(V) for several realizations at <j> = 0.55 and 
F D = 0.025. 



from ten simulations in which the probe particle is se- 
lected at different locations in the sample and driven over 
equal distances. In this work all P(V) curves are normal- 
ized such that the sum of the probabilities equals 1.0. In 
Fig. 3(b) a large peak appears at V = 0.025, correspond- 
ing to time periods when the probe particle moves freely 
without colliding with any background particles, as in- 
dicated by the upper plateaus at V = 0.025 in V(t) in 
Fig. 3(a). Figure 3(b) shows that there is a second peak 
in P(V) centered at V = 0.0125, which is half of the 
value of the free particle velocity. This peak corresponds 
to the time the probe particle spends pushing one addi- 
tional disk, which doubles the drag. The same double 
peak feature in P(V) appears for all <fi < 0.5. In this 
regime, we find that as <j> increases, the relative weight 
of the V = 0.025 peak decreases and is shifted to the 
V = 0.0125 peak. In Fig. 3(c,d) we plot V(t) and P(V) 
for a higher <\> = 0.55. Here the weight at V = 0.025 is 
significantly reduced and the peak probability has shifted 
to V — 0.008. There is still a peak at V = 0.0125 which 
is now higher than than the V = 0.025 peak. 

For <f> = 0.71, the system is dense enough that the 
probe particle spends very little time moving freely out 
of contact with other particles. In Fig. 4(a) we illustrate 
that V{t) for = 0.71 has an intermittent character with 
large velocity spikes. A blow up in Fig. 4(b) of one of the 
velocity spikes from Fig. 4(a) shows that there is some 
asymmetry in the pulses. A slow build up in the velocity 
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FIG. 4: V(t) for (p = 0.71 and F D = 0.025, showing a more 
intermittent motion, (b) A blowup of one of the velocity 
pulses indicating an asymmetry in the response. 




FIG. 5: P(V) for <j> = 0.71 and F D = 0.025, showing a 
peak at low V near V = 0.037 with a long tail extending 
up to V = 0.025. (b) The same P(V) on a log-linear plot 
indicating that the tail has an exponential form. 




time time 



FIG. 6: V(t) for cj> = 0.807 and F D = 0.025, showing increas- 
ing intermittency. (b) A blowup of one of the velocity pulses 
showing the same asymmetry found for <fi — 0.71. 
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FIG. 7: (a) P(V) for <f> = 0.807 at F D = 0.025 showing a 
peak probability at V = 0.0025 with a long tail, (b) The 
same P(V) on a log-linear scale, showing that the tail fits to 
an exponential form. 



is followed by a sharper velocity drop. In Fig. 5(a), P(V) 
for <j) = 0.71 has a maximum near V = 0.00375 with 
a long tail ending at V = 0.025. There are only very 
small peaks near V = 0.0125 and V = 0.025 which are 
remnants of the type of dynamics seen in Fig. 3(a,b). In 
Fig. 5(b) we plot P(V) on a log-linear scale to indicate 
that the tail has an exponential form. There are some 
deviations from a single exponential fit for V > 0.015 due 
to the existence of an upper bound on the velocities which 
generates a peak at V = 0.025. The intermittency and 
asymmetry of the velocity jumps are very similar to the 
observations in two-dimensional drag experiments near 
the jamming transition 23]. In those experiments, the 
stress rather than velocity fluctuations were measured, so 
the asymmetry is reversed. The experiments of Ref. [23j 
also showed that the fluctuations had a maximum with an 
exponential tail feature similar to that seen in Fig. 5(b) 
In Fig. 6(a) we show V{t) for <fi = 0.807, where the in- 
termittency is more pronounced. The maximum value of 
the velocity spikes is always below V = 0.025, indicating 
that the probe particle is always in contact with the sur- 
rounding particles and never moves freely. The increasing 
intermittency is similar in appearance to the stress fluctu- 
ations seen in shear experiments 24]. Figure 6(b) shows 
a blowup of a velocity spike from Fig. 6(a) indicating that 
the same asymmetry observed for <f> = 0.71 persists. In 
Fig. 7(a) we plot P(V) for 4> = 0.807, where the peak 
probability is shifted to a lower velocity V — 0.0025 and 
the velocity tail is more pronounced. Fig. 7(b) shows 



P(V) at cf> = 0.807 on a log-linear plot illustrating that 
the tail has an exponential form. In the dragging experi- 
ments of Geng and Behringer in Ref. [23| , the maximum 
density that was studied was around (f> = 0.77, close to 
the (ft = 0.807 density presented in Fig. 7. In the experi- 
mental work, the avalanche events were found to have an 
exponential distribution very similar to the fluctuations 
we observe in Fig. 7(a). The fact that the slow drag 
experimental system in Ref. 23] could not go to higher 
packing fractions is likely due to frictional forces, which 
are absent in our system. In general we observe the same 
type of exponential decay in the tail of the velocities for 
densities up to within two percent of 4>j. 

In Fig. 8 we plot V(t) for (f> = 0.8414, which is within 
a fraction 0.998 of 4>j. Here the fluctuations show a pro- 
nounced intermittency where the velocity drops almost to 
zero between velocity bursts. In Fig. 9(a) we show P(V) 
for (j> = 0.8414 in log- log format and in Fig. 9(b) we show 
the same data on a log-normal scale. The light line is a a 
power law fit with an exponent of a = —3.0. The power 
law fit is valid for slightly more than a decade, so at this 
density we cannot rule out a stretched exponential form; 
however, a purely exponential fit of the type performed 
at the lower densities of 4> = 0.807 and <\> — 0.71 no longer 
fits the data. 

In Fig. 10(a) we plot P(V) for <j) = 0.823, 0.833, 0.8395, 
and 0.8414 on a log-log scale while in Fig. 10(b) we plot 
the same curves on a log-linear scale to show that the 
velocity distribution tails change from an exponential to 
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FIG. 8: V(t) close to jamming at <j> = 0.8414 and F D = 0.025. 
Here V drops almost to zero between the velocity spikes. 




FIG. 9: P(V) at <f> = 0.8414 and F D = 0.025. (a) Log-log 
scale with a power law fit (light line) with a = —3.0. (b) The 
same as in (a) on a log-normal scale, showing that a power 
law fit (light line) is better than an exponential fit. 



a power law form. The densities (j> = 0.823 and 0.833 
are fitted well by an exponential, but this fit deviates for 
4> = 0.8395 and 0.8414, crossing over to a power law form. 
If point J is a critical point as suggested in recent studies, 
than the various quantities should diverge as a power law 
as 4>j is approached and the dynamical response should 
also reflect the critical properties. In general, critical phe- 
nomena or fluctuations should only be relevant within 1 
percent of the critical point or closer, which is consistent 




FIG. 10: P(V) for 4> = 0.823, 0.833, 0.8395, and 0.8414 [from 
lower left to upper left in panel (a)] at Fd = 0.025. The solid 
line is a power law fit with a — —3.0 and the dotted line is 
an exponential fit. (a) Log- log scale, (b) Log- normal scale. 



with our observations. We note that the power law ex- 
ponent of a — —3.0 used to fit our data close to jamming 
is comparable to results obtained for recent experimen- 
tal studies of granular slip systems where the energy is 
released in avalanche events distributed with power law 
exponents of —2.8 [25]. 

Several factors limit how closely we can approach <j)j. 
One is the fact that we have a finite but small drive on 
the particle of Fd = 0.025. Ideally, a vanishingly small 
drive would be best, but with smaller drives prohibitively 
long simulation times are necessary to obtain adequate 
statistics. Another issue is that for densities close to <fij, 
it possible that the probe particle can move for some dis- 
tance and then jam. This effect was noted in the work of 
O'Hcrn et al. who found that the distribution of possible 
jammed configurations peaks around 4> = 0.844 but that 
there was still some finite width to this peak as long as 
the system was finite [2J . 

Velocity fluctuations with power law distributions have 
also been observed in chute flow systems, where three 
phases are identified as a function of density: a low den- 
sity free-falling region, a small liquid-like region, and a 
high density solid or glassy region [2(| . In the liquid-like 
region, the velocities can be fit to a power law with expo- 
nents of —2.4 to —3.7. Even though the chute flow system 
is three-dimensional, one possible interpretation of our 
two-dimensional results is that the low density regime is 
similar to the free-falling regime, while the glassy or solid 
and the liquidlike regimes could correspond to the region 
just before jamming where we find power law distribu- 
tions. 

It is interesting to ask why the non-exponential fea- 
tures such as those in Fig. 9 were not observed in the 
slow drag experiments of Ref. (23[. In the experiments 
the jamming density fell at a lower (f> than our (j)j , which 
is likely due to frictional effects. Although the exper- 
imental system is jammed, it does not appear to have 
the same critical properties that are associated with the 
frictionless grains at <fij = 0.843. This implies that the 
concept of critical jamming may be limited to certain 
types of models of jamming and that criticality at point 
J may occur only in the frictionless limit. 

In Fig. 11 we plot the trajectories of the particles as the 
probe moves halfway through the sample. In Fig. ll(a,b) 
at (p = 0.129 and 0.323, the main trajectory of the probe 
particle appears as a continuous line. Only small per- 
turbations of the surrounding particles appear near the 
probe particle, showing that the perturbations are all 
highly localized. For <f> = 0.807 in Fig. 11(c), the pertur- 
bation extends much further from the probe particle in 
the transverse direction, with motion occurring for par- 
ticles that are up to nine lattice constants away from the 
driven particle. Here the trajectories of the surround- 
ing particles are very short and do not intertwine. For 
<j> = 0.8414 in Fig. 11(d), the extent of the maximum 
transverse perturbations is only slightly larger than at 
4> = 0.807; however, the perturbations ahead of the probe 
particle are much larger. In addition, the motion of the 
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FIG. 11: Black lines: trajectories followed by the particles 
as the probe particle moves halfway through the system at 
F D = 0.025. (a) (j> = 0.129. (b) <j> = 0.323 (c) <f> = 0.807. (d) 
d> = 0.8414. 



surrounding particles takes two different forms. The first 
type of motion consists of those particles which move in 
straight trajectories, as seen at lower <fi. The second type 
of motion involves particles which move in larger circu- 
lar trajectories that intertwine with the trajectories of 
other particles. These trajectories are associated with 
plastic deformations in which some of the surrounding 
particles move around one another. Similar plastic re- 
arrangements of the surrounding particles caused by the 
motion of a probe particle were observed in a charged 
colloidal media 17]. In Ref. [ljj there was a threshold 
for the motion of the driven particle, unlike the motion 
in Fig. 11(d) which occurs with no threshold. The plas- 
tic deformation trajectories that we observe are similar 
to topological changes or Tl events. We also note that 
Fig. 2(b) indicated that (V) does not develop a nonlin- 
ear dependence on cj> until <fi > 0.825. We only observe 
the intersecting plastic trajectories for <f> > 0.825. It is 
possible that the nonlinearity in (V) may be associated 
with the onset of the plastic motions. 



A. Power Spectra 

Another quantity that was measured in the slow drag 
experiments of Ref. [23| is the power spectrum 
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FIG. 12: The power spectra S(u>) of V(t) at constant Fd = 
0.025 for (a) <j> = 0.484, (b) cj> = 0.71, (c) (j> = 0.823 and (d) 
4> = 0.841. For (a-c) a the high frequencies have a ui~ 2 form 
indicated by a solid line. A fit to a Lorenztian is indicated 
by a dashed line in (b) and a similar fit can be performed for 
(a,c). For <j> — 0.841 in (d), a Lorenztian can no longer be fit 
and a low frequency w _1 feature appears. 



In Ref. [23], the force fluctuations showed an w~ 2 behav- 
ior at high frequencies with a knee or roll off to a flat spec- 
trum at low frequencies. The roll off frequency shifted to 
lower frequencies with lower rotation rates. The same 
type of behavior was also observed for increasing 
densities and no clear signature appeared in the power 
spectrum of the fluctuations just before jamming. In our 
system, the frequency range is cut off by the inverse of 
the time required for the particle to transverse half of the 
system, so for higher densities where the probe particle 
moves more slowly, we can access lower frequencies. 

For <f) < 0.84, we find that the power spectra of the 
velocity fluctuations have the same features as the ex- 
periment in Ref. (23[. In Fig. 12(a,b,c) we plot the power 
spectra S(lu) for <fi = 0.484, 0.71, and 0.823 showing an 
lu~ 2 feature at high frequency and a knee feature at lower 
frequencies that is very similar to those seen in experi- 
ments. The power spectra can be fit to a Lorentzian form 
as indicated in Fig. 12(b), 

a/ \ Squjq 
S ( u} ) <* 7a~, — 2- 



S(u) 



v{ty 



- iujt dt 



Here ujq is the crossover frequency at the knee and So is 
the spectral power along the fiat portion of the spectrum 
below the knee. The presence of a Lorentzian can be in- 
terpreted as the emergence of a characteristic time scale. 
At the lower <j> the time scale could be the average time 
between collisions or an average time between the large 
velocity bursts. In general, for increasing <f> we find that 
the knee feature shifts to lower frequencies. 

In Fig. 12(d) we plot S{w) for <j> = 0.841 within the 
regime where the velocity histograms showed a power 
law form. Here the spectrum no longer has a Lorentzian 
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FIG. 13: (a) P(V) at <f> = 0.8414 for F D = 0.025, 0.125, and 
0.5, from top left to bottom left, (b) P(V) at <j> = 0.8426 for 
F D = 0.06, 0.125, and 0.25, from top left to bottom left. 

form but shows al/u signal at low frequencies followed 
by the persistent 1/lu 2 feature at higher frequencies. The 
appearance of l/u> noise is indicative of the presence of 
scale- free fluctuations, which are consistent with a critical 
phenomenon. We note that if we were able to perform 
very long time traces and access lower frequencies, the 
low frequency 1 /lu noise power would be higher than the 
noise power in the flat portion of the Lorentzian found 
at lower densities, implying that the low frequency noise 
power So may diverge as the jamming transition is ap- 
proached. A divergence in the noise power at generic sec- 
ond order phase transitions has recently been proposed 
[27l | which is consistent with our results. We also note 
that a peak in the noise power and the onset of 1// 
noise were observed two-dimensional systems at disor- 
dering transitions [28| . 

The slow drag experiments of Ref. [23| may not have 
been able to access the true critical properties of point J 
since the fluctuations were examined at a lower 4> than 
the expected <j)j. The drag experiments of Albert et al. 
[2l| also produced Lorentzian type power spectra. We 
note that this was in a three-dimensional system where 
frictional effects may also have been relevant. 

It would be interesting to see whether a crossover from 
Lorentzian to 1// noise occurs experimentally as jam- 
ming is approached. In Ref. [2{|, a probe particle was 
pushed from the bottom of a packing of rodlike particles 
and a transition from l// 2 noise to 1// noise was corre- 
lated with a transition from stick-slip type motion to the 
motion of a solid plug. Although this result is consistent 
with our results, the experimental system is somewhat 
different in that there are no particle rearrangements in 
the solid plug state and the noise is due to friction effects 
with the wall. 



B. Driving Rate Dependence 

We next examine how the velocity fluctuations vary 
with the magnitude of the driving force Fd on the probe 
particle. In Fig. 13 we plot P(V) at <j) = 0.8414 for 
F D = 0.025, 0.125, and 0.5. For F D = 0.125 there is 
a peak in P(V) at V = 0.015 and for F D = 0.5 the 
peak shifts to V — 0.1. In both cases the tail also shifts 
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FIG. 14: (a) Circles: (V) vs \<j>j - <f>\/(f>j for F D = 0.025. 
Straight line: power law fit with a — 2.0. (b) Circles: (V) vs 
\4>J ~ 4>\/^J f° r Fd = 0.125. Straight line: power law fit with 
a = 1.75. 

to higher V. For Fd = 0.125 the velocity tail can still 
be fit to a power law; however, the range of the fit is 
reduced from the Fd = 0.025 case and for Fd = 0.5 the 
range of the fit is even smaller. This suggests that the 
critical fluctuations are most prominent in the limit of 
small driving forces. 

In Fig. 13(b) we plot P(V) on a log-log scale at 
F D = 0.06, 0.125, and 0.25 for a system with = 0.8426 
which is within 0.001 of the jamming density. At this 
density, for Fd < 0.06 it is difficult to obtain enough 
statistics for a histogram due to the very long simulation 
times; also, when the system is this close to the jamming 
density, the probe particle can in some cases jam before 
it moves halfway through the system. At the jamming 
density the probe particle almost aways jams instantly. 
For the lowest drive of Fd = 0.06, Fig. 13(b) shows that 
P{V) has a power law form. The data can be fit with an 
exponent of a — —2.7 over a range of two decades, while 
if we use only data with V > 0.03, a fit with a = —3.1 is 
better. Although we cannot more firmly establish the ex- 
act exponent, Fig. 13(b) indicates that the distributions 
cannot be fit well to a purely exponential form. 

C. Estimation of the Correlation Length Scale 
Exponent v 

In Ref. , a measurement of the number of disks that 
were pushed along with the probe particle was used to 
obtain the exponent of a diverging correlation length, 

Zk(</>j- </>)-". (l) 

The resulting value of v was between 0.6 and 0.7. Other 
estimates of the exponent for the same two-dimensional 
disk system have given v = 0.71 ± 0.12 v = 0.6 [lj|, 
and v = 0.57 [T3|. More recent quasistatic shearing sim- 
ulations [3(3] have found higher estimates of v = 0.8 and 
v = 1.0; it is also argued that the exponent is sensitive 
to the exact value of <j>j, which may explain why the 
exponents in Refs. [!, Qj(| are somewhat smaller. Refer- 
ences @,[l(| both used about 1200 particles, with Ref. [|[ 
reporting <pj = 0.8395 for a driven probe particle and 
Ref. [l(| giving 4>j = 0.8415 for a shearing simulation. 
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In Ref. [121] . where a higher value of cj>j = 0.8433 was 
obtained, about 2500 particles were used. In the present 
work we have a comparable number of particles N « 2600 
and we find 4>j = 0.843 in agreement with Ref. [Hj]. More 
recent shearing simulations with larger systems than in 
1C | have found a higher value of 4>J = 0.843 and v = 0.8 
30l | so we believe a value of (j>j = 0.843 is reasonably 



accurate. 

In Fig. 2(b), the velocity (V) drops to zero in a nonlin- 
ear fashion as <fr increases toward <j)j and the concavity 
in the velocity also changes very near jamming. We can 
estimate v based on the scaling behavior of (V) as the 
jamming transition is approached. At <j> = 0, the probe 
particle moves with a free velocity of Vo- At nonzero 
</>, the velocity of the probe particle decreases accord- 
ing to y k Vo/N m , where N m is the number of other 
disks that are being pushed along by the probe parti- 
cle. If the pushed particles are within a disklikc jammed 
area of radius £ that surrounds the probe particle, then 
N m = 4>j^ 2 ■ If £ oc (</>j — then the velocity should 

scale as V oc Vo(<j) — 4>j) 2v ■ From the velocity versus 
cj> curves we can fit V oc (</> — 4>j) a , giving v — a/2. 
In Fig. 14(a) we plot V vs \<f)j - cj)\/<j>j for F D = 0.025 
along with a fit of a = 2.0, which gives v = 1.0. Figure 
14(b) shows the same plot for Fd = 0.125 where we ob- 
tain a = 1.75 and v — 0.875. These values are within 
the range of the estimates from the quasistatic shearing 
simulations of Ref. We note that other theoretical 

studies have found v = 0.5 [H, IS HI- Although our scal- 
ing range is too small to give a definitive value for v ? our 
results along with the other simulations 0, OH H, G3 
are consistent with v > 0.5. It is possible that there 
is more than one diverging length scale near jamming. 
From the images of the particle displacements near jam- 
ming in Fig. 11, there appear to be two different modes 
of displacements of the surrounding particles: the small 
displacements that extend in front of the driven particle 
and the plastic deformations that occur further from the 
probe particle. It may be that there is a different length 
scale associated with each of these two types of motion. 



D. Velocity Force Curves 




FIG. 15: (V) vs Fd (dark curves) and d{V)/dF D vs F D (light 
curves). The velocity-force curves are linear for (a) <j> = 0.483 
and (b) cj> — 0.807 and are nonlinear for (c) <j> = 0.8414 and 
(d) 4> = 0.8426. 
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FIG. 16: (V) vs F D on a log-log scale for </> = 0.483 (tri- 
angles), 4> — 0.807 (diamonds), (j> — 0.8414 (squares), and 
4> = 0.8426 (circles). The straight lines are power law fits with 
= 1.0 for 4> = 0.483 and 0.807, f3 = 1.24 for <j> = 0.8414, and 
f3 = 1.47 for 4> = 0.8426. 



We next examine the average velocity (V) versus Fjj 
for varied (p. At <fi — in the absence of other disks, the 
overdamped dynamics we employ cause the velocity of 
the probe particle to vary linearly with the applied force, 
V oc Fd- For tj> well below <f>j, collisions with other 
particles reduce the probe particle velocity and increase 
the effective viscosity; however, the velocity-force curves 
remain linear as shown in Fig. 15(a,b) for 4> = 0.483 and 
= 0.807. The plots of d{V)/dF D in Fig. 15(a,b) have 
little curvature, which is consistent with a linear velocity- 
force relationship. Very close to jamming the velocity- 
force curves become nonlinear as seen in Fig. 15(c,d) for 
(j) = 0.8414 and <f> = 0.8426, where the d(V)/dF D curves 
show growing deviation from linearity. 



In Fig. 16 we plot (V) vs Fd. on a log-log scale where 
the solid lines are power law fits with the form V = F^. 
For (f) = 0.483 and (j) = 0.807 the fits are linear with 
0=1, while for = 0.8414 and = 0.8426, the fits 
give = 1.24 and = 1.47, respectively. In general, 
exponents are not expected to change continuously near 
a critical point so the exponent at <j> = 0.8414 probably 
originates from a smeared combination of linear and non- 
linear regions. The exponent of = 1.47 at <f> = 0.8426 
should be more accurate and is also in agreement with a 
number of previous studies of single driven probe parti- 
cles. Simulations of a driven particle moving through a 
disordered assembly of charged colloids gave an exponent 
of = 1.5 [IliGil- In both the colloidal system and the 
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disk system, there are plastic rearrangements of the par- 
ticles around the probe particle; however, there are sev- 
eral differences between the two systems, including the 
fact that the colloidal particle is effectively jammed at 
low drives and has a nonzero critical threshold force for 
motion. Additionally, in the charged colloid system there 
is no well defined jamming density <j)j due to the longer 
range soft Yukawa interaction potential. Experiments 
where a single probe particle is driven through a three- 
dimensional colloidal system as the density is increased 
towards the glass transition gave an exponent of (3 = 1.5 
at densities below the glass transition, which can be in- 
terpreted as below the jamming transition [lj|. The clos- 
est experimental realization of our system is the work in 
Ref. [161 ] where a probe particle is driven through a two- 
dimensional assembly of grains. One difference between 
the experiments in Ref. [16| and our simulations is that 
the experiment contains an additional vibration to induce 
a fluidized state. In the experiment, linear velocity-force 
curves appear at low densities, while closer to jamming 
there is a nonlinear region at low drives with a crossover 
to a linear region at high drives. This crossover shifts 
to higher values of Fjj as the jamming transition is ap- 
proached. In the nonlinear regime, Ref. [16[ used a fit 
to V oc exp(Fjj) rather than a power law; however, the 
range of the data is sufficiently small that it does not rule 
out a power law fit. Additionally, it is possible that the 
applied vibrations induce activation dynamics that are 
not present in our work. Despite these differences, the 
onset of nonlinear velocity- force curves and increased in- 
termittency close to jamming found in Ref. [l6j are con- 
sistent with our results. 

We note that a power law fit to the velocity force curves 
has been associated with criticality at the onset of motion 
in other two-dimensional systems. In particular, exten- 
sive studies have been performed on plastic depinning of 
particle assemblies driven over rough landscapes in sys- 
tems ranging from colloidal particles moving over disor- 
dered substrates [3l|, [HJ , vortices in type-II superconduc- 
tors [33j , and generalized models of particle motion over 
quenched disorder [U- In these systems, the velocity 
response can be fit to the form 



V=(F D - F c f. (2) 



Here Fc is the threshold force that must be applied for 
motion to occur. For (f> < (j)j, we find Fq = 0; however, 
velocity-force curves with nonlinear scaling can still oc- 
cur even in the absence of a threshold. In the collective 
depinning systems, exponents of (3 = 1.5 to (3 = 2.25 are 
generally observed. Our results showing a crossover from 
linear to nonlinear velocity-force curves near jamming are 
a further indication that point J exhibits critical proper- 
ties and that this critically also appears in transport. 




FIG. 17: (a) (V) vs F D below jamming at <j> = 0.8426 (circles) 
and above jamming at (f> = 0.8439 (squares), (f> = 0.8446 
(diamonds), and (f> = 0.846 (triangles). Above jamming, a 
threshold force F c is needed to depin the probe particle, (b) 
The threshold force F c vs <j> with point J denoting where the 
system becomes jammed. 



IV. THRESHOLD AND YIELDING ABOVE 
JAMMING 

We next show that for <j> > <j>j, there is a finite thresh- 
old force required to move the probe particle through 
the background of disks. In the jammed phase, it is 
known that the onset of motion or yielding consists of 
a coexistence between an unjammed fluid phase and a 
jammed solid phase, and that in the presence of a shear, 
the motion may be confined within a certain region due 
to shear banding [35j. In Fig. 17(a) we plot (V) vs Fd for 
<t> = 0.8426, 0.8439, 0.8446, and 0.846. Here the probe 
is immobile or pinned at low drives for <f> > <fij as indi- 
cated by the fact that (V) = 0.0 until the external force 
Fd reaches a critical force F c for motion. At densities 
just below jamming, Fig. 16(b) showed that the velocity- 
force curves have positive curvature with (3 > 1. Above 
jamming, the velocity-force curves have a different form 
and are more consistent with either a linear dependence 
j3 = 1 or with (3 < 1.0. Although we do not have the ac- 
curacy in this regime to determine the exact form of the 
velocity-force curves above jamming, our results suggest 
that the onset of motion is different for <fi > (fij than for 

The threshold force F c increases with increasing cj>. 
In Fig. 17(b) we plot F c vs (j> showing the onset of the 
threshold at point J and the subsequent rapid increase 
of F c with increasing cf>. Well above <fij, F c increases 
roughly linearly with increasing <j). The appearance of 
a threshold and the general shape of the F c curve is in 
agreement with recent three-dimensional colloid simula- 
tions in which a threshold for motion appears above a 
certain packing density [20]. We note that in our sys- 
tem, the disks interact with a harmonic potential and are 
soft enough that the probe particle can squeeze between 
neighboring disks above jamming. A similar yielding ap- 
pears in shearing simulations performed above jamming 
ficl [l2| . For real granular systems with frictional effects 
and much stiffer interactions between particles, it would 
be difficult to move the probe particle through the sys- 
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FIG. 18: Particle trajectories during the time period required 
for the probe particle to move halfway through the system at 
a drive just above the depinning threshold F c . (a) cj> — 0.8439. 
(b) cj> = 0.846. 



tem above jamming. For colloidal particles with much 
softer interactions, a threshold force can be extracted as 
has been demonstrated in experiment [il^ ]. 

In Fig. 18(a,b) we plot the particle trajectories at a 
drive just above the depinning threshold for <p = 0.8439 
and (j) = 0.846, where plastic deformations occur within 
a localized region. The trajectories are qualitatively sim- 
ilar to the case 4> = 0.8414 shown in Fig. 11(d) where 
the motion consisted of a combination of particles close 
to the driven particle being pushed linearly over a small 
distance while particles further from the driven particle 
undergo intertwining plastic motion. One noticeable dif- 
ference is that the spatial extent of the perturbations 
above jamming is less than that seen below jamming at 
<\> = 0.8414. Additionally, the spatial extent of the per- 
turbations above jamming does not appear to rapidly 
grow in size as 4> increases above <j)j. The fact that the 
motion is more localized above jamming suggests that in 
this regime the system is only plastically distorted near 
the probe while the remaining portion of the system acts 
rigidly elastic. We call this a local shear banding effect, 
and it suggests that the onset of motion or yielding above 
jamming is consistent with the behavior of a disordered 
elastic solid. A similar conclusion that the state above 
jamming is a disordered elastic solid was drawn from two- 
dimensional shearing simulations [i~2| . 

In Fig. 19 we plot P(V) for <j> = 0.846 at F D = 1.725, 
1.85, 4.0, and 6.0. At F D = 1.725 and F D = 1.85, 
P{V) has a broad form with non-exponential features. 
Attempts to perform a single power law fit for the higher 
values of V in this regime give only a limited range with 
an exponent of —3.0. For the higher drives Fd = 4.0 and 
Fq = 6.0, there is still a tail at high drives; however, a 
power law fit cannot be performed. 

In Fig. 20(a) we show the power spectrum S(ui) of V(t) 
at (j> = 0.846 for Fd = 1.725, just above depinning. The 
solid line is a power law fit to S(ui) tx uj~°. We obtain 
a = 0.8 over a wide range of low frequencies, with a 
small region of lo~ 2 dependence at high frequencies. The 
smallness of the lo~ 2 regime is a reversal from what we 




FIG. 19: Velocity histograms P(V) above jamming at cj> = 
0.846 for F D = 1.725, 1.85, 4.0, and 6.0, from upper left to 
lower left. 
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FIG. 20: The power spectra S(u>) of the velocity time traces 
V{t) at 4> = 0.846. (a) F D = 1.725, just above depinning. 
The solid line is a power law fit with a = —0.8 (b) Fd = 6.0. 
A flat spectral signature appears at low frequencies. 



observe below jamming, where the oj~ 2 frequency depen- 
dence dominates over a large range until just below the 
jamming transition. In Fig. 20(b) we plot S(ui) for the 
same density (f> — 0.846 but at Fd = 6.0, well above the 
depinning threshold. In this case, at low frequencies a 
white spectrum with a = appears. 



V. SUMMARY 

In summary, we have numerically investigated the ve- 
locity fluctuations of a driven probe particle as it moves 
through a bidisperse assembly of disks in two dimensions. 
The particular model we study is composed of harmon- 
ically interacting disks with a 1 : 1.4 radius ratio. We 
find a jamming density of cj>j = 0.843, in agreement with 
recent numerical work on the same system (l2l [Til [l6| . 
For a constant driving force on the probe particle, we find 
that the average probe velocity goes to zero continuously 
as (j> approaches <fij from below and that the velocity fluc- 
tuations develop increasingly intermittent features close 
to <j)j. The velocity distribution functions have long tails 
which fall off exponentially, similar to the force fluctu- 
ations observed in the slow drag experiments of Geng 
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and Behringer (23[. When (f> is within one percent of 4>j, 
the velocity distributions can no longer be fit to a purely 
exponential form and have a power law shape suggest- 
ing that the system is exhibiting fluctuations of a criti- 
cal type. The extent of the perturbations caused by the 
driven grain in the surrounding packing also grows as 
jamming as approached and the perturbations take two 
forms. The first is small linear displacements of less than 
a lattice constant, while the second is large winding dis- 
placements which are associated with plastic rearrange- 
ments. The plastic rearrangements occur only very near 
point J where the power law velocity distributions are 
observed. The power spectrum of the fluctuations be- 
low jamming has a Lorentzian shape with a knee at a 
crossover frequency and u>~ 2 behavior at higher frequen- 
cies, similar to the spectral signatures found in slow drag 
experiments [23j . Just below the jamming transition, the 
power spectrum changes shape and develops a l/u> low 
frequency behavior. The velocity-force curves are lin- 
ear below jamming, but become nonlinear with a power 
law form just at the jamming transition. The appear- 
ance of nonexponential velocity distributions, the onset 
of 1 /ui noise and the development of power law velocity- 



force curves are all consistent with jamming having the 
signatures of a critical phenomenon. From the velocity 
versus 4> curves, we obtain an estimate for the correla- 
tion length exponent v = 0.8 — 1.0, in agreement with 
recent studies [12l |30| . For <j) > (f> j we find that a critical 
threshold force F c is required to set the probe particle 
in motion, in agreement with recent numerical studies 
of driven probe particles in dense colloidal assemblies at 
densities above jamming [l8j]. The threshold force in- 
creases with increasing cf> and the motion consists of lo- 
calized plastic distortions near the moving particle which 
we call a local shear banding effect. The velocity fluctu- 
ations at drives just above the threshold force have a 
nonexponential form and the power spectrum exhibits a 
1/uj shape. For higher drives, the velocities are exponen- 
tially distributed, suggesting that the onset of yielding 
by the probe particle for cf> > <f>j has similar properties 
to yielding at point J. 
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useful comments. This work was carried out under the 
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